Chapter Nine 


Frequency Domain Controller Design 


9.1 Introduction 

Frequency domain techniques together with the root locus method have been very 
popular classical methods for both analysis and design of control systems. As 
they are still used extensively in industry for solving controller design problems 
for many industrial processes and systems, they have to be included in academic 
curricula and modern textbooks on control systems. This chapter is organized 
as follows. 

In Section 9.2 we study the open- and closed-loop frequency transfer func¬ 
tions and identify the frequency response parameters such as system frequency 
bandwidth, peak resonance, and resonant frequency. 

In Section 9.3 we show how to read the phase and gain stability margins, and 
the values of the steady state error parameters K p , K v , K a from the corresponding 
frequency diagrams, known as Bode diagrams. Bode diagrams represent the 
magnitude and phase plots of the open-loop transfer function with respect to the 
angular frequency to. They can be obtained either analytically or experimentally. 
In this chapter we present only the analytical study of Bode diagrams, though 
it should be mentioned that Bode diagrams can be obtained from experimental 
measurements performed on a real physical system driven by a sinusoidal input 
with a broad range of frequencies. 

The controller design technique based on Bode diagrams is considered in 
Section 9.4. It is shown how to use the phase-lag, phase-lead, and phase-lag- 
lead controllers such that the compensated systems have the desired phase and 
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gain stability margins and steady state errors. It is also possible in some cases 
to improve the transient response since the damping ratio is proportional to the 
phase margin, and the response rise time is inversely proportional to the system 
bandwidth. 

Section 9.5 contains a case study for a ship positioning control system, and 
Section 9.7 represents a laboratory experiment. In Section 9.6 we comment on 
discrete-time controller design. 

Chapter Objectives 

The main objective of this chapter is to show how to use Bode diagrams as 
a tool for controller design such that compensated systems have, first of all, the 
desired phase and gain stability margins and steady state errors. Controllers based 
on Bode diagrams can also be used to improve some of the transient response 
parameters, but their design is far more complicated and far less accurate than 
the design of the corresponding controllers based on the root locus method. 


9.2 Frequency Response Characteristics 

The open- and closed-loop system transfer functions are defined in Chapter 2. 
For a feedback control system these transfer functions arc respectively given by 

Open-loop: G(s)H(s), Closed-loop : ^ + g(s)H(s) = (9 I) 

The frequency transfer functions arc defined for sinusoidal inputs having all 
possible frequencies t o E [0,+oc). They are obtained from (9.1) by simply 
setting s = ju>, that is 

Open-loop: G{ju)H(ju% Closed-loop: ^ + G(ju)H(ju) = 

(9.2) 

Using the frequency transfer functions in the system analysis gives complete 
information about the system's steady state behavior, but not about the system’s 
transient response. That is why controller design techniques based on frequency 
transfer functions improve primarily the frequency domain specifications such as 
phase and gain relative stability margins and steady state errors. However, some 
frequency domain specifications, to be defined soon, can be related to certain 
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time domain specifications. For example, the wider the system bandwidth, the 
faster system response, which implies a shorter response rise time. 

Typical diagrams for the magnitude and phase of the open-loop frequency 
transfer function are presented in Figure 9.1. From this figure one is able to read 
directly the phase and gain stability margins and the corresponding phase and 
gain crossover frequencies. In Section 9.3, where we present the Bode diagrams, 
which also represent the magnitude and phase plots of the open-loop frequency 
transfer function with respect to frequency, with magnitude being calculated in 
decibels (dB), we will show how to read the values for K v , K v . K a . 
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Figure 9.1: Magnitude (a) and phase (b) of the open-loop transfer function 

In addition to the phase and gain margins, obtained from the open-loop fre¬ 
quency transfer function, other frequency response parameters can be obtained 
from the frequency plot of the magnitude of the closed-loop frequency trans¬ 
fer function. This plot is given in Figure 9.2. The main closed-loop frequency 
response parameters are: system bandwidth, peak resonance, and resonant fre¬ 
quency. They are formally defined below. 










384 


FREQUENCY DOMAIN CONTROLLER DESIGN 


System Bandwidth: This represents the frequency range in which the 
magnitude of the closed-loop frequency transfer function drops no more than 
3 dB (decibels) from its zero-frequency value. The system bandwidth can be 
obtained from the next equality, which indicates the attenuation of 3 dB, as 

\M(juBw)\ = -^=\M(0)\ =>• u>BW (9.3) 

It happens to be computationally very involved to solve equation (9.3) for 
higher-order systems, and hence the system bandwidth is mostly determined 
experimentally. For second-order systems the frequency bandwidth can be found 
analytically (see Problem 9.2). 

Peak Resonance: This is obtained by finding the maximum of the function 
|with respect to frequency lo. It is interesting to point out that the 
systems having large maximum overshoot have also large peak resonance. This 
is analytically justified for a second-order system in Problem 9.1. 

Resonant Frequency: This is the frequency at which the peak resonance 
occurs. It can be obtained from 

^|M(jw)| = 0 u r 



Figure 9.2: Magnitude of the closed-loop transfer function 
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9.3 Bode Diagrams 

Bode diagrams are the main tool for frequency domain controller design, a 
topic that will be presented in detail in Section 9.4. In this section we show 
how to plot these diagrams and how to read from them certain control system 
characteristics, such as the phase and gain stability margins and the constants 
required to determine steady state errors. 

Bode diagrams represent the frequency plots of the magnitude and phase of 
the open-loop frequency transfer function G(jut)H (jut). The magnitude is plotted 
in dB (decibels) on the logo; scale. In general, the open-loop frequency transfer 
function contains elementary frequency transfer functions representing a constant 
term (static gain) and dynamic elements like system real poles and zeros and 
complex conjugate poles and zeros. We first study independently the magnitude 
and frequency plots of each of these elementary frequency transfer functions. 
Since the open-loop frequency transfer function G(jut)H(jut) is given in terms 
of products and ratios of elementary transfer functions, it is easy to see that 
the phase ofG(jut)H (jut) is obtained by summing and subtracting phases of the 
elementary transfer functions. Also, by expressing the magnitude of the open- 
loop transfer function in decibels, the magnitude \G(jut)H(jut)\ dB is obtained 
by adding the magnitudes of the elementary frequency transfer functions. For 
example 


\G(jut)H(jut)\ dB - 20 log 10 


K(jut + z x )(jut + z>) 
{jut)(jut + p 2 )(jut + p 3 ) 


- 20 log 10 \K\ + 20 log 10 | jut + %! | + 20 log 10 | jut + 2 2 


and 


+ 20 log 10 


1 

jut 


+ 20 log 10 


1 

jut + p-2 


+ 20log 10 


1 

jut + P3 


arg {G(jut)H(jut)} - arg {K} + arg {jut + z x j + arg {jut + z 2 } 

- arg {jut} - arg {jut + p 2 } - arg {jut + p 3 } 

In the following we show how to draw Bode diagrams for elementary 
frequency transfer functions. 
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Constant Term: Since 


K d B = 20 log 10 K 


positive number K > 1 
negative number K < 1 


r o°, k > o 

\ —180°, K < 0 


(9.4) 


the magnitude and phase of this element are easily drawn and arc presented in 
Figure 9.3. 





0 arg{/f>0) 


► logco 


- 180 ° 


arg{/T<0) 


Figure 9.3: Magnitude and phase diagrams for a constant 
Pure Integrator: The transfer function of a pure integrator, given by 

G{M = — (9.5) 

has the following magnitude and phase 

|G(jw )| dB = 20 login - = —20log 10 lj, arg G(ju) = -90° (9.6) 

CO 

It can be observed that the phase for a pure integrator is constant, whereas the 
magnitude is represented by a straight line intersecting the frequency axis at 
to = 1 and having the slope of — 20 dB/decade. Both diagrams arc represented 
in Figure 9.4. Thus, a pure integrator introduces a phase shift of — 90° and a 
gain attenuation of —20 dB/decade. 
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Figure 9.4: Magnitude and phase diagrams for a pure integrator 


Pure Differentiator: The transfer function of a pure differentiator is given 
by 

G{ju) = ju (9.7) 

Its magnitude and phase arc easily obtained as 

\G{jw)\ dB ~ 20log 10 w, arg G(ju) = 90° (9.8) 


The corresponding frequency diagrams are presented in Figure 9.5. It can be 
concluded that a pure differentiator introduces a positive phase shift of 90° and 
an amplification of 20 dB/decade. 

Real Pole: The transfer function of a real pole, given by 


G CH 


p 

P + ju 


1 

i+7f 


(9.9) 


has the following magnitude and phase 
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Figure 9.5: Magnitude and phase diagrams for a pure differentiator 


l G OOU = - 201 °gio 


1 + 



1/2 


? 


arg G(ju>) = — tan 1 



(9.10) 


The phase diagram for a real pole can be plotted directly from (9.10). It can 
be seen that for large values of lo, lo > p, the phase contribution is —90°. For lo 
small, lo <C p. the phase is close to zero, and for lo = p the phase contribution is 
—45°. This information is sufficient to sketch arg G(jio) as given in Figure 9.6. 

For the magnitude, we see from (9.10) that for small lo the magnitude is 
very close to zero. For large values of lo we can neglect 1 compared to c o/p 
so that we have a si mi lar result as for a pure integrator, i.e. we obtain an 
attenuation of 20dB/decade. For small and large frequencies we have straight- 
line approximations. These straight lines intersect at lo = p, which is also known 
as a comer frequency. The actual magnitude curve is below the straightline 
approximations. It has the biggest deviation from the asymptotes at the corner 
frequency (see Figure 9.6). 
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Figure 9.6: Magnitude and phase diagrams for a real pole 

Real Zero: The transfer function of an element representing a real zero is 
given by 


G(ju) = 7O + jw) = 1 + j(^p) 


(9.11) 


Its magnitude and phase arc 


\G(ju)\ dB = 201og 10 



arg G(ju >) = tan 1 



(9.12) 


Using analysis similar to that performed for a real pole, we can conclude that 
for small frequencies an asymptote for the magnitude is equal to zero and for 
large frequencies the magnitude asymptote has a slope of 20 dB/decade and 
intersects the real axis at lo = z (the corner frequency). The phase diagram for 
small frequencies also has an asymptote equal to zero and for large frequencies 
an asymptote of 90°. The magnitude and phase Bode diagrams for a real-zero 
element are represented in Figure 9.7. 
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Figure 9.7: Magnitude and phase diagrams for a real zero 

Complex Conjugate Poles: The transfer function of an element representing 
a pair of complex conjugate poles is in fact a transfer function of a second-order 
system, which has the form 

GCH = - 


1 


OO 2 + 2(l - 

The magnitude and phase of this second-order system arc given by 


\°{j^)\dB = —20 log 10 


2£cu 


UJr 


+ (1-^ 

loz 


2\2 


1/2 


axg G(ju>) = — tan 


-l 




For large values of lo the coiTesponding approximations of (9.14) arc 

,2 ' 


|G(jw)| 


d.B 


LO 


LOt 


-20log 10 ( — = —40 log 10 I — 


LO 


LO r 


(9.13) 


(9.14) 
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arg {<?(ju?)} « — tan 1 — — tan 1 (0 ) = —180° 

At low frequencies the approximations can be obtained directly from (9.13), that 

18 , 2 

G(jco) « = 1 =* I G{jL0)\ dB = 0, arg {G(ju>)} = 0° 

n 

Thus, the coiTesponding asymptotes for small and large frequencies are, respec¬ 
tively, zero and —40dB/decade (with the corner frequency at lo = u> n ) for the 
magnitude, and zero and —180° for the phase. At the corner frequency u> n the 
phase is equal to —90°. The coiTesponding Bode diagrams are represented in 
Figure 9.8. Note that the actual plot in the neighborhood of the corner frequency 
depends on the values of the damping ratio £. Several curves are shown for 
0.1 < ( < 1. It can be seen from Figure 9.8 that the smaller (, the higher peak 
of the magnitude plot. 



Figure 9.8: Magnitude and phase diagrams for complex conjugate poles 
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Complex Conjugate Zeros: An element that has complex conjugate zeros 
can be represented in the form 


G{ju) 


1 + 2 O' 




(9.15) 


so that the corresponding Bode diagrams will be the mirror images of the Bode 
diagrams obtained for the complex conjugate poles represented by (9.13). In the 
case of complex conjugate zeros, the asymptotes for small frequencies are equal to 
zero for both the magnitude and phase plots; for high frequencies the magnitude 
asymptote has a slope of 40 dB /decade and stalls at the corner frequency of 
lo — Lo n , and the phase plot asymptote is 180°. 


9.3.1 Phase and Gain Stability Margins from Bode Diagrams 

It has been already indicated in Figure 9.1 how to read the phase and gain stability 
margins from the frequency magnitude and phase plots of the open-loop feedback 
transfer function. In the case of Bode diagrams the magnitude plot is expressed 
in dB (decibels) so that the gain crossover frequency is obtained at the point 
of intersection of the Bode magnitude plot and the frequency axis. Bearing in 
mind the definition of the phase and gain stability margins given in (4.54) and 

(4.55) , and the corresponding phase and gain crossover frequencies defined in 

(4.56) and (4.57), it is easy to conclude that these margins can be found from 
Bode diagrams as indicated in Figure 9.9. 

Example 9.1: In this example we use MATLAB to plot Bode diagrams for 
the following open-loop frequency transfer function 




1 

T-1 

+ 

3 

jv{ju> + 2) 

( j u> ) + 2(ytc) + 2 


Bode diagrams are obtained by using the MATLAB function bode (num, den) . 
The phase and gain stability margins and the phase and gain crossover frequencies 
can be obtained by using [Gm, Pm, wep, weg] =marg±n (num, den) . Note 
that the open-loop frequency transfer function has to be specified in terms of 
polynomials num (numerator) and den (denominator). The MATLAB function 
conv helps to multiply polynomials as explained below in the program written 
to plot Bode diagrams and find the phase and gain margins for Example 9.1. 
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Figure 9.9: Gain and phase margins and Bode diagrams 


num=[1 1] ; 

dl= [ 1 0]; 

d2= [ 1 2]; 

d3=[1 2 2 ]; 

denl=conv(dl,d2); 

den=conv(deni, d3) ; 

bode(num,den); 

[Gm,Pm,wcp,wcg]=marg±n(num,den); 

The corresponding Bode diagrams are presented in Figure 9.10. The phase 
and gain stability margins and the corresponding crossover frequencies arc ob¬ 
tained as 

Gm = 8.9443dB, Pm = 82.2462°, lo cv = 1.7989rad/s 9 u> cg = 


0.2558 rad/s 
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Example 9.2: Consider the following open-loop frequency transfer function 


G{juj)H{juj) 


10 (ju; + 4) 

{ju + l)(i w + + 3) 


The corresponding Bode diagrams obtained by following the same MATLAB 
instructions as in Example 9.1 arc given in Figure 9.11. The phase and gain 
margins and the phase and gain crossover frequencies are obtained as 


Gm = oc, Pm = 43.1488°, io cp = oo, u> cg = 3.0576rad/s 


■o 


9.3.2 Steady State Errors and Bode Diagrams 

Steady state errors can be indirectly determined from Bode diagrams by reading 
the values for constants K p , K v , K a from them. Knowing these constants, the 
corresponding errors arc easily found by using formulas (6.30), (6.32), and (6.34). 
The steady state errors and corresponding constants K p . K v . K a are first of all 
determined by the system type, which represents the multiplicity of the pole at 
the origin of the open-loop feedback transfer function, in general, represented by 
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Figure 9.11: Bode diagrams for Example 9.2 


G(Jw)H(ju>) 

This can be rewritten as 


K(jio + Zi)(ju} + z 2 ) ■■■ 

(ju> + Pi)(ju> +p 2 ) ■■■ 


G{ju)H{jJ) 


™2"-W(i + £)(i + £) ••• 


^(i + f)(i + ff) ••• 
M'(i + ff)(i + £) ••• 


(9.16) 


(9.17) 


where 

Kz 1 z 2 ••• 

K b = - 

VlP2 ••• 

is known as Bode’s gain, and r is the type of feedback control system. 


(9.18) 
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For control systems of type r — 0, the position constant according to formula 
(6.31) is obtained from (9.17) as 


A- B (l+f)(l+g) ■■■ 

w"(i + f)(i + £) lj “=° 


(9.19) 


It follows from (9.17)—(9.19) that the corresponding magnitude Bode diagram of 
type zero control systems for small values of lo is flat (has a slope of 0 dB) and the 
value of 20 log Kb — 20 log K v . This is graphically represented in Figure 9.12. 


\G(j(n)H(j(o)\ 



Figure 9.12: Magnitude Bode diagram of type 
zero control systems at small frequencies 


For control systems of type r = 1. the open-loop frequency transfer function 
is approximated at low frequencies by 


*«( i +£)( i +e0- k b 

(M'^ + frXi + g)-" 


(9.20) 


It follows that the corresponding magnitude Bode diagram of type one control 
systems for small values of iv has a slope of —20 dB/decade and the values of 


20 log 


Kb 

jw 


— 20 log \Kb \ — 20 log |cu 


(9.21) 


From (9.20) and (6.33) it is easy to conclude that for type one control systems 
the velocity constant is K v — Kb- Using this fact and the frequency plot of 
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(9.21), we conclude that K v is equal to the frequency at which the line (9.21) 
intersects the frequency axis, that is 

0 = 20 log K B - 20 log |w* | K B = w* = IC V (9.22) 

This is graphically represented in Figure 9.13. 



Figure 9.13: Magnitude Bode diagram of type 
one control systems at small frequencies 


Note that if K v — to* > 1, the corresponding frequency i a* is obtained at the 
point where the extended initial curve, which has a slope of — 20 dB/decade, 
intersects the frequency axis (see Figure 9.13b). 

Similarly, for type two control systems, r — 2, we have at low frequencies 


*«(*+ £)(■•+£) - Kb 

UK (i + ff) (i + £)••• " 


(9.23) 


which indicates an initial slope of —40 dB/decade and a frequency approximation 
of 


20 log 


K b 

U w f 


— 20log \K b \ - 20log cu 2 = 20log \K B \ — 40logcu 


(9.24) 


From (9.23) and (6.35) it is easy to conclude that for type two control systems 
the acceleration constant is K a — I{ B . From the frequency plot of the straight 
line (9.24), it follows that K B — (cu**) 2 , where uj** represents the intersection of 
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the initial magnitude Bode plot with the frequency axis as represented in Figure 
9.14. 



Figure 9.14: Magnitude Bode diagram of type 
two control systems at small frequencies 


It can be seen from Figures 9.12-9.14 that by increasing the values for 
the magnitude Bode diagrams at low frequencies (i.e. by increasing Kb), the 
constants K v . K „, and K a arc increased. According to the formulas for steady 
state errors, given in (6.30), (6.32), and (6.34) as 


1 


1 + K v ’ 


1 

K~, 


2 

^SSp a ra bol ic 7V" 

-ft a 


we conclude that in this case the steady state errors arc decreased. Thus, the 
bigger Kb, the smaller the steady state errors. 

Example 9.3: Consider Bode diagrams obtained in Examples 9.1 and 9.2. 
The Bode diagram in Figure 9.10 has an initial slope of —20 dB/decade which 
intersects the frequency axis at roughly w* = 0.2 rad/s. Thus, we have for the 
Bode diagram in Figure 9.10 


K p = oo, K v S3 0.2, K a = 0 


Using the exact formula for K v , given by (6.33), we get 

(« + l> 


K = lim < s 


s—>0 [ s(s + 2)(s 2 + 2s + 2) 


= 0.25 
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In Figure 9.11 the initial slope is 0 dB, and hence we have from this diagram 
20 log IFj, « 15 =>■ K p « 5.62, K v = 0, K a = 0 
Using the exact formula for K p as given by (6.31) produces 

K p = lim { -- ^ --) = 6.67 

s—*o \ (s + l)(s + 2)(s + 3) / 

Note that the accurate results about steady state error constants arc obtained easily 
by using the corresponding formulas; hence the Bode diagrams are used only for 
quick and rough estimates of these constants. 


o 


9.4 Compensator Design Using Bode Diagrams 

In this section we show how to use Bode diagrams in order to design controllers 
such that the closed-loop system has the desired specifications. Three main types 
of controllers—phase-lead, phase-lag, and phase-lag-lead controllers—have been 
introduced in the time domain in Chapter 8. Here we give their interpretation 
in the frequency domain. 

We present the design procedure for the general controllers mentioned above. 
Similar and simpler procedures can be developed for PD, PI, and PID controllers. 
After mastering the design with phase-lead, phase-lag, and phase-lag-lead con¬ 
trollers, students will be able to propose their own algorithms for PD, PI, and 
PID controllers. 

Controller design techniques in the frequency domain will be governed by the 
following facts: 

(a) Steady state errors are improved by increasing Bode’s gain Kg. 

(b) System stability is improved by increasing phase and gain margins. 

(c) Overshoot is reduced by increasing the phase stability margin. 

(d) Rise time is reduced by increasing the system’s bandwidth. 

However, very often it is not possible to satisfy all of these requirements at 
the same time, and control engineers have to compromise between several 
contradicting requirements. 
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The first two items, (a) and (b), have been already clarified. In order to 
justify item (c), we consider the open-loop transfer function of a second-order 
system given by 


y 2 y J + 20„) 

whose gain crossover frequency can be easily found from 

,2 


\G(ju cg )H(jio cg )\ = 


LIZ 


LOy/uj 2 + 4 ( 2 U 2 

leading to 

Wog=U%\f\TrfW^e 

The phase of (9.25) at the gain crossover frequency is 


= 1 


arg {G(ju} cg )H(ju cg )} = -90° - tan 


-l W c g 


2(U> n 


so that the corresponding phase margin becomes 

2 C 


Pm = tan 


^^TT4e--2C 2 


= Pm (() 


(9.25) 

(9.26) 

(9.27) 

(9.28) 

(9.29) 


Plotting the function Tm((), it can be shown that it is a monotonically increasing 
function with respect to £; we therefore conclude that the higher phase margin, 
the larger the damping ratio, which implies the smaller the overshoot. 

Item (d) cannot be analytically justified since we do not have an analytical 
expression for the response rise time. However, it is very well known from 
undergraduate courses on linear systems and signals that rapidly changing signals 
have a wide bandwidth. Thus, systems that are able to accommodate fast signals 
must have a wide bandwidth. 

In the remainder of this section, we first present standard controllers (phase- 
lag, phase-lead, and phase-lag-lead) in the frequency domain, and then show how 
to use these in order to achieve the desired system specifications. Each design 
technique will be given in an algorithmic form, and each will be demonstrated 
by an example. 
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9.4.1 Phase-Lag Controller 

The transfer function of a phase-lag controller is given by 


Glagiju) = — 


Pi \ Zl + ju 


Zi J Pi + 3^ 


i + JiT 
1 + J’ft' 


Zi > pi 


(9.30) 


The corresponding magnitude and phase frequency diagrams for a phase-lag 
controller are presented in Figure 9.15. 




Figure 9.15: Magnitude approximation and exact phase of a phase-lag controller 

Note that for the magnitude diagram it is sufficient to use only the straightline 
approximations for a complete understanding of the role of this controller. In 
general, straightline approximations can be used for almost all Bode magnitude 
diagrams in controller design problems. However, phase Bode diagrams are very 
sensitive to changes in frequency in the neighborhood of the corner frequencies, 
and so should be drawn as accurately as possible. 

Due to attenuation of the phase-lag controller at high frequencies, the 
frequency bandwidth of the compensated system (controller and system in series) 
is reduced. Thus, the phase-lag controllers are used in order to decrease the 
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system bandwidth (to slow down the system response). In addition, they can be 
used to improve the stability margins (phase and gain) while keeping the steady 
state errors constant. 

Expressions for u> max and o,, iax of a phase-lag controller will be derived in 
the next subsection in the context of the study of a phase-lead controller. As a 
matter of fact, both types of controllers have the same expressions for these two 
important design quantities. 


9.4.2 Phase-Lead Controller 


The transfer function of a phase-lead controller is 


Gl e ad(j( J j') 


(P_2_\ Z-2 + ju _ 1 + 

\Z2jP2+ju t+jg’ 


P2 > Z2 


(9.31) 


and the corresponding magnitude and phase Bode diagrams are shown in Figure 
9.16. 



Figure 9.16: Magnitude approximation and exact phase of a phase-lead controller 

Due to phase-lead controller (compensator) amplification at higher frequencies, 
it increases the bandwidth of the compensated system. The phase-lead controllers 
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are used to improve the gain and phase stability margins and to increase the system 
bandwidth (decrease the system response rise time). 

It follows from (9.31) that the phase of a phase-lead controller is given by 


arg {G lead (ju>)} = tan 




so that 

arg {Glead^ju)} = 0 u> max = s Jz 2 p -2 

Assume that 


P-2 = a > 1 

Substituting Lo max in (9.32) implies 



tail (prnax 


(2 — 1 
2 y/a 


(9.32) 


(9.33) 


(9.34) 


(9.35) 


It is left as an exercise for students to give detailed derivations of formula 
(9.35)—see Problem 9.3. 

It is easy to find, from (9.35), that the value for parameter a in terms of 
4>max is given by 


1 “F sill (f>max 

1 sill tfimax 


(9.36) 


Note that the same formulas for u> m ax, (9.33), and the parameter a, (9.36), 
hold for a phase-lag controller with pi , z\ replacing p >, z- 2 and with p\ = az \, 
a < 1. 


9.4.3 Phase-Lag-Lead Controller 

The phase-lag-lead controller has the features of both phase-lag and phase-lead 
controllers and can be used to improve both the transient response and steady 
state errors. However, its design is more complicated than the design of either 
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phase-lag or phase-lead controllers. The frequency transfer function of the phase- 
lag-lead controller is given by 


G c (jio) 


(j iu + £l)(. Vf_ + ^ 2 ) 
U U + Pl){ju + P 2 ) 


Z\ Z'2 

m ’ 2 ( 1 + ^)( 1 + J s ) 


( 1 + -' k )( 1 + -'£) 


1)^1 = PlP2, P-2 > Z-2 > tl > pi 


(9.37) 


The Bode diagrams of this controller arc shown in Figure 9.17. 



Figure 9.17: Bode diagrams of a phase-lag-lead controller 


9.4.4 Compensator Design with Phase-Lead Controller 

The following algorithm can be used to design a controller (compensator) with 
a phase-lead network. 
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Algorithm 9.1: 


1. Determine the value of the Bode gain Kb given by (9.18) as 

Kz 1 z 2 ■ ■ ■ 

K b = - 

PlP2 ■■■ 

such that the steady state error requirement is satisfied. 

2. Find the phase and gain margins of the original system with Kb determined 
in step 1. 

3. Find the phase difference, A <f>, between the actual and desired phase margins 
and take o max to be 5°-10° greater than this difference. Only in rare cases 
should this be greater than 10°. This is due to the fact that we have to give 
an estimate of a new gain crossover frequency, which can not be determined 
very accurately (see step 5). 


4. Calculate the value for parameter a from formula (9.36), i.e. by using 

1 T sin (f> ma x . 
a = i- ^~Z - > 1 

1 Sill (f>max 


5. Estimate a value for a compensator’s pole such that uj max is roughly located 
at the new gain crossover frequency, to max ~ uj cgn ew ■ As a rule of thumb, 
add the gain of AG = 20log(a)[dB] at high frequencies to the uncompen¬ 
sated system and estimate the intersection of the magnitude diagram with 
the frequency axis, say uq. The new gain crossover frequency is somewhere 
in between the old ui cg and uq. Some authors (Kuo, 1991) suggest fixing 
the new gain crossover frequency at the point where the magnitude Bode 
diagram has the value of —0.5AG[dB]. Using the value for parameter a 
obtained in step 4 find the value for the compensator pole from (9.34) as 
— p c — —Umaxs/a and the value for compensator’s zero as — z c — —p c /a. 
Note that one can also guess a value for p c and then evaluate z c and u> max . 
The phase-lead compensator now can be represented by 


G c {s) 


as + Pc 

S+Pc 


6. Draw the Bode diagram of the given system with controller and check the 
values for the gain and phase margins. If they arc satisfactory, the controller 
design is done, otherwise repeat steps 1-5. 
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The next example illustrates controller design using Algorithm 9.1. 
Example 9.4: Consider the following open-loop frequency transfer function 


G(ju>)H(ju>) 


_ K(ju + 6) _ 

{ju + l)(i w + 2 -){ju> + 3) 


Step 1. Let the design requirements be set such that the steady state error due to 
a unit step is less than 2% and the phase margin is at least 45°. Since 

1 1 K X 6 

e ss = -= -, Afi= - = A 

1 + K v 1 + AV 1x2x3 

we conclude that K > 50 will satisfy the steady state error requirement of being 
less than 2%. We know from the root locus technique that high static gains 
can damage system stability, and so for the rest of this design problem we take 
K = 50. 

Step 2. We draw Bode diagrams of the uncompensated system with the Bode gain 
obtained in step 1 and determine the phase and gain margins and the crossover 
frequencies. This can be done by using the following sequence of MATLAB 
functions. 

[den]=±nput('enter denominator'); 

% for this example [den]=[l 6 11 6]; 

[num]=input('enter numerator'); 

% for this example [num]=[50 300]; 

[Gm,Pm,wcp,wcg]=marg±n(num,den); 
bode(num,den) 

The corresponding Bode diagrams arc presented in Figure 9.18a. The phase 
and gain margins arc obtained as Gm = oo, Pm — 5.59° and the crossover 
frequencies are u) cg = 7.5423rad/s, u> cp = oo. 

Step 3. Since the desired phase is well above the actual one, the phase-lead 
controller must make up for 45° — 5.59° = 39.41°. We add 10°, for the 
reason explained in step 3 of Algorithm 9.1, so that 4> max = 49.41°. The above 
operations can be achieved by using MATLAB as follows 
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Figure 9.18: Bode diagrams for the original system 
(a) and compensated system (b) of Example 9.4 


% estimate Phimax with Pmd = desired phase margin; 
Pmd=input('enter desired value for phase margin'); 
Phimax=Pmd-Pm+l0; 

% converts Phimax into radians; 

Ph±rad=(Phimax/180)*pi; 

Step 4. Here we evaluate the parameter a according to the formula (9.36) and 
get a — 7.3144. This can be done in MATLAB by 

a=(l + sin(Phirad))/(1-sin(Phirad) ) ; 

Step 5. In order to obtain an estimate for the new gain crossover frequency 
we first find the controller amplification at high frequencies, which is equal to 
20 log (a) = 17.2836 dB = A Gas- The magnitude Bode diagram increased by 
A GdB at high frequencies intersects the frequency axis at cu,„ :ax . ~ 10.5 rad/s. 
We guess (estimate) the value for p c as p c — 25, which is roughly equal to 
^maxVti- By using p c — 25 and forming the corresponding compensator, we get 
for the compensated system Prric — 48.2891° at u> cgnew — 13.8519rad/s, which 
is satisfactory. This step can be performed by MATLAB as follows. 
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% Find amplification at high frequencies, DG; 
DG=20*logl0(a) ; 

% estimate value for pole —pc from Step 5; 
pc=input('enter estimated value for pole pc'); 

% form compensator's numerator; 
nc=[a pc]; 

% form compensator's denominator; 
dc=[1 pc]; 

% find the compensated system transfer function; 
numc=conv(num,nc); 
denc=conv(den,dc) ; 

[Gmc,Pmc,wcp,wcg]=margin(numc,dene) ; 
bode(numc,dene) 

The phase-lead compensator obtained is given by 

7.3144s + 25 as + p -2 
s + 25 s + p -2 

Step 6. The Bode diagrams of the compensated control system arc presented in 
Figure 9.18b. Both requirements arc satisfied, and therefore the controller design 
procedure is successfully completed. 

It is interesting to compare the transient response characteristics of the 
compensated and uncompensated systems. This cannot be easily done analytically 
since the orders of both systems arc greater than two, but it can be simply 
performed by using MATLAB. Note that num, den, numc, dene represent, 
respectively, the numerators and denominators of the open-loop transfer functions 
of the original and compensated systems. In order to find the corresponding 
closed-loop transfer functions, we use the MATLAB function cloop, that is 

[enum, eden]=cloop(num,den,-1); 

% —1 indicates a negative unit feedback 
[enume,cdenc]=cloop(numc, dene, -1) ; 

The closed-loop step responses are obtained by 

[y,x]=step(enum, eden) ; 

[yc,xc]=step(enume, cdenc) ; 

and arc represented in Figure 9.19. It can be seen from this figure that both 
the maximum percent overshoot and the settling time arc drastically reduced. In 


G e {s) 
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addition, the rise time of the compensated system is shortened since the phase- 
lead controller increases the frequency bandwidth of the system. 



Figure 9.19: Step responses for the original (a) and compensated (b) systems 

The highly oscillatory behavior of the step response of the original system 
can be fully understood from its root locus, which is given in Figure 9.20. It can 
be seen that the real part of a pair of complex conjugate poles is very small for 
almost all values of the static gain, which causes high-frequency oscillations and 
very slow convergence to the response steady state value. 

The closed-loop eigenvalues of the original and compensated systems, ob¬ 
tained by MATLAB functions roots (cden) and roots (cdenc) , are given 
by 

Ai = -5.3288, A 2 , 3 = -0.3356 ± j7.5704 
Ai,2c = -9.6736 ± jf 13.3989* A 3c = -8.2628, A 4c = -3.3900 

which indicates a big difference in the real parts of complex conjugate poles for 
the original and compensated systems. 


o 
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Figure 9.20: Root locus of the original system 

9.4.5 Compensator Design with Phase-Lag Controller 

Compensator design using phase-lag controllers is based on the compensator’s 
attenuation at high frequencies, which causes a shift of the gain crossover 
frequency to the lower frequency region where the phase margin is high. The 
phase-lag compensator can be designed by the following algorithm. 

Algorithm 9.2: 

1. Determine the value of the Bode gain K B that satisfies the steady state error 
requirement. 

2. Find on the phase Bode plot the frequency which has the phase margin 
equal to the desired phase margin increased by 5° to 10°. This frequency 
represents the new gain crossover frequency, i o cgntw . 

3. Read the required attenuation at the new gain crossover frequency, i.e. 
\AG(jLO cg new)\ dB , and find the parameter a from 

— 20 log = —20 log (a) = |A G(j 

U C gnew)\ c iB 
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which implies 


a 




1 

\AG(joj 

cgnew ) | 


Note that 


|A G(jio 

cgnew ) | 


\K 11 cgnew H" %1 \ \ j^ cgnew + *2 
| j^cgnew "1“ Pi \ \ cgnew + P 2 I • 


4. Place the controller zero one decade to the left of the new gain crossover 
frequency, that is 

_ ^cgnew 

" c ~ 10 


Find the pole location from p c 
compensator has the form 

G c (s) 


az c = au) cgnew /W. The required 

as + p c 
s+ Pc 


5. Redraw the Bode diagram of the given system with the controller and check 
the values for the gain and phase margins. If they are satisfactory, the 
controller design is done, otherwise repeat steps 1-5. 


Example 9.5: Consider a control system represented by 

G(s) = - -^- 

v ' s(s + 2)(s + 30) 

Design a phase-lag compensator such that the following specifications are 
met: e SSrarnF < 0.05, Pm > 45°. The minimum value for the static 
gain that produces the required steady state error is equal to K = 1200. 
The original system with this static gain has phase and gain margins 
given by Pm = 6.6449°, Gm = 4.0824dB and crossover frequencies of 
u; cg = 6.1031 rad/s, u) cp = 7.746 rad/s. 

The new gain crossover frequency can be estimated as c o cgn ew = 1.4 rad/s 
since for that frequency the phase margin of the original system is approximately 
50°. At c Ocgnew = 1.4 rad/s the required gain attenuation is obtained by 
MATLAB as 
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wcgnew=l.4; 
dl=1200; 
gl=abs ( j *1.4); 
g2=abs(j *1.4 + 2); 
g3=abs(j *1.4 + 30) ; 
dG=dl/(gl*g2*g3); 

which produces AG(jl. 1) = 11.6906 and a — l/|AG(jl.4)| = 0.0855. The 
compensator’s pole and zero are obtained as — z c — —LV cgnew /10 = —0.14 and 
— p c — —auj C g new / 10 — —0.0120 (see step 4 of Algorithm 9.2). The transfer 
function of the phase-lag compensator is 

, 0.0855s + 0.0120 

G c (s) = - 

v ; s +0.0120 

The Bode diagrams of the original and compensated systems arc given in Figure 
9.21. 



KT' 10° 10 1 10 2 

Frequency (rad/sec) 


Figure 9.21: Bode diagrams for the original system 
(a) and compensated system (b) of Example 9.5 

The new phase and gain margins and the actual crossover frequencies arc 
I'rnc = 47.03°, Gmc = 24.82 dB, u cgnew = 1.405 rad/s, cc cpJieu , = 7.477 rad/s 
and so the design requirements arc satisfied. The step responses of the original 












FREQUENCY DOMAIN CONTROLLER DESIGN 


413 


and compensated systems are presented in Figure 9.22. 



Figure 9.22: Step responses for the original system 
(a) and compensated system (b) of Example 9.5 


It can be seen from this figure that the overshoot is reduced from roughly 0.83 to 
0.3. In addition, it can be observed that the settling time is also reduced. Note 
that the phase-lag controller reduces the system bandwidth ( u cgnew < u cg ) so 
that the rise time of the compensated system is increased. 


o 

9.4.6 Compensator Design with Phase-Lag-Lead Controller 

Compensator design using a phase-lag-lead controller can be performed according 
to the algorithm given below, in which we first form a phase-lead compensator 
and then a phase-lag compensator. Finally, we connect them together in series. 
Note that several different algorithms for the phase-lag-lead controller design 
can be found in the control literature. 

Algorithm 9.3: 

1. Set a value for the static gain Kb such that the steady state error requirement 
is satisfied. 

2. Draw Bode diagrams with Kb obtained in step 1 and find the corresponding 
phase and gain margins. 
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3. Find the difference between the actual and desired phase margins, 
A <f> — Pmd — Pm, and take (f) max to be a little bit greater than Ao. Cal¬ 
culate the parameter a 2 of a phase-lead controller by using formula (9.36), 
that is 

1 T sin (fimax 

a 2 = T-- 

1 Sill (pmax 

4. Locate the new gain crossover frequency at the point where 

20 log \G(jui cgnew )\ - -101oga 2 (9.38) 


5. Compute the values for the phase-lead compensator’s pole and zero from 


Pc2 — ^ cgnew V«2, Zc2=Pc2/a-2 (9.39) 

6. Select the phase-lag compensator’s zero and pole according to 

z c i - 0.1z c2 , p c 1 - z c i/a 2 (9.40) 


7. Form the transfer function of the phase-lag-lead compensator as 


G c (s) — Glag(s) X Gl ea d(s ) 


Zcl 


Pel 


8 + Z c 2 
s +Pc 2 


8. Plot Bode diagrams of the compensated system and check whether the design 
specifications are met. If not, repeat some of the steps of the proposed 
algorithm—in most cases go back to steps 3 or 4. 

The phase-lead part of this compensator helps to increase the phase margin 
(increases the damping ratio, which reduces the maximum percent overshoot and 
settling time) and broaden the system’s bandwidth (reduces the rise time). The 
phase-lag part, on the other hand, helps to improve the steady state errors. 

Example 9.6: Consider a control system that has the open-loop transfer 
function 

G( s ) - _ _ + 10 ^ _ 

W (s 2 + 2s + 2)(s + 20) 

For this system we design a phase-lag-lead controller by following Algorithm 
9.3 such that the compensated system has a steady state error of less than 4% 
and a phase margin greater than 50°. In the first step, we choose a value for 
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the static gain K that produces the desired steady state error. It is easy to check 
that K — 100 e ss — 3.85%, and therefore in the following we stick with this 
value for the static gain. Bode diagrams of the original system with K — 100 
are presented in Figure 9.23. 


Gm=lnf dB, (w= NaN) Pm=31.61 deg. (w=7.668) 



Frequency (rad/sec) 


Figure 9.23: Bode diagrams of the original system 

It can be seen from these diagrams—and with help of MATLAB determined 
accurately—that the phase and gain margins and the corresponding crossover 
frequencies are given by Pm — 31.61°, Gm — oo, and u cg — 7.668rad/s, 
uj cp = oo. According to step 3 of Algorithm 9.3, a controller has to introduce 
a phase lead of 18.39°. We take 4> max — 25° and find the required parameter 
a- 2 — 2.4639. Taking u cgnew — 20 rad/s in step 4 and completing the design 
steps 5-8 we find that Pm — 39.94°, which is not satisfactory. We go back to 
step 3 and take o, nax = 30° = 0.5236 rad, which implies u> — 3. 

Step 4 of Algorithm 9.3 can be executed efficiently by MATLAB by per¬ 
forming the following search. Since —10log 3 = —10.9861 dB we search the 
magnitude diagram for the frequency where the attenuation is approximately 
equal to —11 dB. We start search at u> = 20 rad/s since at that point, according 
to Figure 9.23, the attenuation is obviously smaller than —11 dB. The following 
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MATLAB program is used to find the new gain crossover frequency, i.e. to 
solve approximately equation (9.38) 

w=2 0; 

while 20*logl0(100*abs(j*w+10)/ 

abs(((j *w)~ 2 + 2 * j *w+2)*(j *w+2 0)))<-l1; 

w=w-l; 

end 

This program produces c o cgnew — 10 rad/s. In steps 5 and 6 the phase-lag-lead 
controller zeros and poles arc obtained as —p c 2 = —17.3205, — z C 2 — —5.7735 
for the phase-lead part and —p c \ — —0.1925, —z c i — —0.5774 for the phase-lag 
part; hence the phase-lag-lead controller has the form 

n ^ _ s + 0.5774 ^ s + 17.3205 
_ s + 0.1925 X s + 5.7735 

The Bode diagrams of the compensated system arc given in Figure 9.24. 

Gm=lnf dB, (w= NaN) Pm=56.34 deg. (w=4.738) 

50 1 . . . . . . . . 
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Figure 9.24 Bode diagrams of the compensated system 


It can be seen that the phase margin obtained of 56.34° meets the design require¬ 
ment and that the actual gain crossover frequency, 4.738 rad/s, is considerably 
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smaller than the one predicted. This contributes to the generally accepted inac¬ 
curacy of frequency methods for controller design based on Bode diagrams. 

The step responses of the original and compensated systems are compared in 
Figure 9.25. The transient response of the compensated system is improved since 
the maximum percent overshoot is considerably reduced. However, the system 
rise time is increased due to the fact that the system bandwidth is shortened 
(^cgnew — 4.738rad/s < u cg — 7.668rad/s). 



Figure 9.25 Step responses of the original (a) and compensated (b) systems 


o 


9.5 MATLAB Case Study 

Consider the problem of finding a controller for the ship positioning control 
system given in Problem 7.5. The goal is to increase stability phase margin 
above 30°. The problem matrices arc given by 


'-0.0546 

0 

0.5435' 


0 

1 

0 

0 

, B = 

0 

0 

0 

-1.55 _ 


1,55 _ 


A = 


, C = [0 1 0], D = 0 
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The transfer function of the ship positioning system is obtained by the 
MATLAB instruction [num, den] =ss2tf (A, B, C, D) and is given by 


w s(s + 1.55)(s + 0.0546) 

The phase and gain stability margins of this system arc Pm = —19.94° and 
Gm = —15.86 dB, with the crossover frequencies uj cji = 0.2909rad/s and 
uj C g = 0.7025 rad/s (see the Bode diagrams in Figure 9.26). From known values 
for the phase and gain margins, we can conclude that this system has very poor 
stability properties. 


Gm=-15.86dB, (w= 0.2909) Pm=-19.94 deg. (w=0.7025) 



Frequency (rad/sec) 


Figure 9.26: Bode diagrams of a ship positioning control system 

Since the phase margin is well below the desired one, we need a controller 
which will make up for almost a 50° increase in phase. In general, it is hard to 
stabilize systems that have large negative phase and gain stability margins. In the 
following we will design phase-lead, phase-lag, and phase-lag-lead controllers to 
solve this problem and compare the results obtained. 

Phase-Lead Controller: By using Algorithm 9.1 with <p max — 50° + 10° = 
60° we get a phase margin of only 23.536°, which is not satisfactory. It is 
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necessary to make up for o n ,, lx — 50° + 27° = 87°. In the latter case the 
compensator has the transfer function 


G c {$) = 76.31 


a + 0.2038 
a + 15.55 


Figure 9.27 shows Bode diagrams of both the original (a) and compensated (b) 
systems. 




Figure 9.27: Bode diagrams for a ship positioning system: 
(a) original system, (b) phase-lead compensated system 


The gain and phase stability margins of the compensated system arc found from 
the above Bode diagrams as Grac — 15.1603 dB, Pmc — 30.0538°, and the 
crossover frequencies are t j cpc — 1.7618rad/s, u cgc — 4.6419rad/s. The step 
response of the compensated system exhibits an overshoot of 45.47% (see Figure 
9.28). 

Phase-Lag-Lead Controller: By using Algorithm 9.3 we find the compen¬ 
sator transfer function as 


_ , N „,s + 0.1599 _ „ s + 0.016 

GJa) = 7.524—- x 0.1329 - 

° v > oriQ 


s + 1.203 


a + 0.002125 
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Time (secs) 

Figure 9.28: Step response of the compensated 
system with a phase-lead controller 

The Bode diagrams of the original and compensated systems are shown in 
Figure 9.29. 



10- 1 10 ° 10 ' 10 2 
Frequency (rad/sec) 


Figure 9.29: Bode diagrams for a ship positioning control system: 
(a) original system, (b) phase-lag-lead compensated system 
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The phase and gain margins of the compensated system arc given by 
Pmc = 39.6694°, Gmc = 14 dB and the crossover frequencies arc 

l o cgc = 0.4332 rad/s, uj cpc = 1.2401 rad/s. 

From the step response of the compensated system (see Figure 9.30), we can 
observe that this compensated system has a smaller overshoot and a larger rise 
time than the system compensated only by the phase-lead controller. 



Figure 9.30: Step response of the compensated 
system with a phase-lag-lead controller 

Phase-Lag Controller: If we choose a new gain crossover frequency at 
cocgnew = 0.03 rad/s, the phase margin at that point will clearly be above 50°. 
Proceeding with a phase-lag compensator design, according to Algorithm 9.2, 
we get |AG(j0.03)| = 290.7390 and a = 0.034, which implies 2 C = 0.003 and 
p c = 1.0319 X 10 5 . Using the corresponding phase-lag compensator produces 
very good stability margins for the compensated system, i.e. Gm = 32.91 dB 
and Pm — 54.33°. The maximum percent overshoot obtained is much better 
than with the previously used compensators and is equal to MPOS = 18%. 
However, the closed-loop step response reveals that the obtained system is too 
sluggish since the response peak time is t p = 95.2381 s (note that in the previous 
two cases the peak time is only a few seconds). 
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One may try to get better agreement by designing a phase-lag compensator, 
which will reduce the phase margin of the compensated system to just above 
30°. In order to do this we write a MATLAB program, which searches the phase 
Bode diagram and finds the frequency corresponding to the prespecified value 
of the phase. That frequency is used as a new gain crossover frequency. Let 
Pm = 35° = 0.6109 rad. The MATLAB program is 


w=0.1; 
while pi+ 

angle(1/((j *w)*(j *w+l.55)*(j*w+0.0546) ))<0.610 9; 

w=w-0. 01; 

end 

dG=0.8424*abs(1/((j*w)*(j*w+1.55)*(j*w+0.0546) )) ; 

This program produces i o cgnew = 0.07rad/s and |AG(j0.07)| = 87.3677. From 
step 4 of Algorithm 9.2 we obtain the phase-lag controller of the form 


G c (s) 


s + 0.009 
s + 0.000081 


x 0.0114 


The Bode diagrams of the compensated system are given in Figure 9.31. 


Gm=23.48 dB, (w= 0.2651) Pm=31.44 deg. (w=0.06043) 




Figure 9.31: Bode diagram of the phase-lag compensated system 
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It can be seen that the phase and gain margins are satisfactory and given by 
Pm = 31.44° and Gm = 23.48 dB. The actual gain crossover frequencies arc 
W cgnew = 0.06043rad/s and u cpnew = 0.2651 rad/s. 

The closed-loop step response of the phase-lag compensated system, given 
in Figure 9.32, shows that the peak time is reduced to t p = 50.15 s—which 
is still fairly big—and that the maximum percent overshoot is increased to 
MPOS = 45.82%, which is comparable to the phase-lead and phase-lag-lead 
compensation. 



Figure 9.32: Step response of the phase-lag compensated system 

Comparing all three controllers and their performances, we can conclude that, 
for this particular - problem, the phase-lag compensation produces the worst result, 
and therefore either the phase-lead or phase-lag-lead controller should be used. 


9.6 Comments on Discrete-Time Controller Design 

Bode diagrams were originally introduced for studying continuous-time systems 
(Bode, 1940). However, discrete-time systems can be studied using the same 
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diagrams. The bilinear transformation, already used in this book in the stability 
study of discrete-time systems, that has the form 


s 


z - 1 

ITT’ 


^ = 


1 + 3 

1 — S 


(9.41) 


maps the imaginary axis from the s -plane into the unit circle in the 2 -plane 
and vice versa. Since on the unit circle 2 = e^ zT , it is easy to establish the 
relationship between angular frequencies in the s and 2 domains. It is left as an 
exercise for students to show that 


w, = tan 


lo 7 T 


lo z = — tan u> s 


(9.42) 


The above transformation allows one to map the discrete-time open-loop 
transfer function into the continuous-time open-loop transfer function, that is 


G{z)\ z= ^_ =G(a) (9.43) 

and to perform controller design in the continuous-time domain. The results 
obtained have to be mapped back into the discrete-time domain by using (9.42), 
that is 


G c (s) ls= ^ = G c (z) (9.44) 

Note that several bilinear transformations, which arc just scaled versions of (9.41), 
can be found in the control literature. For more details the reader is referred, 
for example, to Franklin et al. (1990), DiStefano etal. (1990), and Phillips and 
Nagle (1995). MATLAB discrete-time controller design problems can be found 
in Shahian and Hassul (1993). 

9.7 MATLAB Laboratory Experiment 


Part 1. Consider the closed-loop system represented in Figure 9.33. This 
system has a transport lag-element, e~ Ts , which represents a time delay of T 
time units. The transport lag-element can be approximated for small values of 
time delay T by 


0 ) 


-Ts 


1 


P) 


.-Ts 


1 + Ts 7 


1 - Ts 
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Figure 9.33: Block diagram of a control system 


Using the following values for the time delay T — 0, 0.01, 0.1, design phase- 
lead, phase-lag, and phase-lag-lead compensators to meet the following closed- 
loop design requirements: steady state error e SSrarnr < 0.02 and Pm > 50°. 
Consider both approximations (i) and (ii) and compare the results obtained. 

Part 2. Draw the exact phase Bode diagrams of the systems compensated 
in Part 1 including the exact contribution from the time delay element. Note that 
the time delay element does not affect the magnitude Bode diagram, but modifies 
the phase Bode diagrams by a factor of /_e~^ T . Compare the approximated and 
exact phase Bode diagrams. Draw conclusions about the impact of time delay 
elements on phase and gain stability margins. 

Part 3. Consider the controller design problem for a system represented by 
its open-loop transfer function 

G{ s ) - _ + 6 ) _ 

1 J (> + 10)(> 2 + 2 s + 2 ) 

This system has been studied in Example 8.9 and in the MATLAB laboratory 
experiment for the root locus controller design in Section 8.8. 

(a) Design a phase-lag-lead controller using Bode diagrams such that the com¬ 
pensated system has the same specifications as those in Section 8.8, i.e. the 
steady state error is less than 1% and the phase margin is such that t s < 2 s 
and MPOS < 10%. Note that the maximum percent overshoot and settling 
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time are inversely proportional to the phase margin. Experiment with sev¬ 
eral values for the phase margin and take the one that satisfies both transient 
response requirements. 

(b) Compare the results obtained with those from Section 8.8 and comment 
on the differences between root locus and Bode diagram phase-lag-lead 
controller design. Which one is easier to design? Which one is more 
accurate? 
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9.9 Problems 


9.1 


Show that for a second-order closed-loop system 

,2 


M(ju) = — 


UJZ 


( j W ) 2 + 2^4^ (j'w) + ^ 

the resonant frequency is given by 


2 

n 


UJ r — sj 1 
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and the peak resonance is 


p 2Cx/r^e 

9.2 Show that the frequency bandwidth for a second-order closed-loop system 
given in Problem 9.1 is 

lobw = U n \j(l- 2C 2 ) + V^C 4 - 4C 2 

Since the 5%-settling time is given by formula (6.20) as 


"5 ' ~ A 

( ,u n 

conclude that the settling time is inversely proportional to the system 
bandwidth, in other words, the wider the system bandwidth, the shorter 
the settling time. 

9.3 Derive formula (9.35) for the maximum phase of a phase-lead controller. 

9.4 Using MATLAB, draw Bode diagrams for a magnetic tape control system 
considered in Problem 5.12. Matrices A and B are given in Problem 5.12. 
The output matrices are 

C = [1 0 1 0], D = 0 

Find the phase and gain stability margins for this system. 

9.5 Based on Algorithms 9.1-9.3, propose algorithms for controller design with 

(a) a PD controller; 

(b) a PI controller; 

(c) a PID controller. 

9.6 Solve the controller design problem defined in Example 9.5 by using both 
phase-lead and phase-lag-lead controllers. 

9.7 Design a phase-lag-lead network for the system 

GO) = . A 

+ 3) 

such that Pm > 45° and e SSrarnp < 0.03. 
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9.8 The block diagram of a servo control system is shown in Figure 9.34. 



Figure 9.34: Block diagram for Problem 9.8 

Design phase-lag, phase-lead, and phase-lag-lead controllers such that the 
phase margin is greater than 60°. 

9.9 A unit feedback system has the transfer function 

rf v. = 40iil 

s(s + 4)0 + 10) 

(a) Construct Bode diagrams for K = 1. 

(b) Using the MATLAB function margin, determine the phase and 
gain margins of the system. 

(c) Determine the value for the static gain K such that the system has 
a gain margin of Gm = 10 dB. Find the corresponding steady state 
errors. 

(d) From Bode diagrams determine the value for the static gain K such 
that the phase margin is 45°. Determine the damping ratio and the 
natural frequency for the obtained value of K. 

9.10 A zero type plant has the transfer function 

GO) = ---- 

0 + 5)0 + 20 ) 

(a) Determine the damping ratio and the natural frequency of the corre¬ 
sponding closed-loop system. 
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(b) Find the steady state errors and the system overshoot. 

(c) Determine the phase and gain margins. 

(d) Design a controller such that the compensated system has a phase 
margin of at least 50° and a steady state error less than 0.02. 

9.11 A unit feedback system with transport delay has the transfer function 

_ 10A'e-°- 5s 

" (7 + 4)(s + 10) 

Approximate time delay as e“ Ts ~ 1/(1 + As). 

(a) Plot on the same figure the Bode diagrams for the system with and 
without transport delay for K = 1. Comment on system stability. 

(b) Plot the unit step response of the system for both cases and compare 
the results (steady state errors, transient response parameters). 

(c) Determine the value of the static gain K that gives a steady state unit 
step error of 0.02 for both the approximated system and the original 
system without time delay. 

(d) Design phase-lag, phase-lead, and phase-lag-lead controllers to 
achieve a steady state error of e ss < 0.05 and a phase margin 
Pm > 50° for the approximated system. 

(e) Draw the exact phase Bode diagrams of the compensated systems, 
including the time delay, and check the values obtained for the phase 
margins. 

9.12 Determine a passive cascade compensator for a unit feedback system 

4 

G(s) = —-- 

s(s 2 + 4s + 8) 

such that the compensated system has a phase margin of 45° and a steady 
state ramp error of less than 2%. 

9.13 Derive formulas (9.42). 

9.14 Solve the controller design problem defined in Example 9.4 by using a 
phase-lag-lead controller. What is the advantage of the controller’s phase- 
lag part? 

9.15 Consider a control system that has the open-loop transfer function 

Q( s ) - _ A (s + 20) _ 

{ } (s + l)(s + 3)(s + 10) 
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(a) Use MATLAB to design any controller by using Bode diagrams such 
that the compensated system has the best possible transient response 
and steady state specifications. 

(b) Solve the same problem using the root locus technique for controller 
design. 

(c) Compare the obtained results and comment on the simplicity (or com¬ 
plexity) of the root locus and Bode diagram methods for controller 
design. 

9.16 Repeat Problem 9.15 for the open-loop control system defined in Problem 
9.9. 

9.17 In order to relate the phase margin (gain margin) and the real parts of the 
system eigenvalues as quantities for relative stability measure, perform the 
following experiment. Consider the system 

G(s) = _ — ^ + _ 

{ ’ {s + 10)(s 2 + 2s + 2) 

(a) Vary the value for static gain K from zero to 100 in increments of 
10, and for each value of K find the phase and gain margins. Plot 
both phase and gain margins as functions of K. 

(b) For each value of K find the closed-loop eigenvalues and plot the 
magnitude of the eigenvalue real parts with respect to K. 

(c) Compare diagrams obtained in (a) and (b) and draw the correspond¬ 
ing conclusion. 

9.18 Repeat Problem 9.17 for the open-loop control system defined in Problem 
9.9. 

9.19 Consider a unit feedback control system that has the open-loop transfer 
function 

= C2 + sy- Ts 

1 (1 + a)(a + 10) 

Note that the term e~ Ts represents a time delay. 

(a) Assuming that the time delay is negligible, draw the corresponding 
Bode diagrams and determine the phase and gain stability margins. 

(b) Since the time delay affects only the phase diagram, draw the 
corrected phase Bode diagram for 7 = 0.1. Determine the phase 
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and gain stability margins and compare them to the corresponding 
quantities found in (a). 

(c) Repeat part (b) for T — 0.2 and T — 0.5. Comment on the impact 
of the time delay element on the phase and gain stability margins. 



